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ABSTRACT 

The character of the first galaxies at redshifts z ^ 10 strongly depends on their level of 
pre-enrichment, which is in turn determined by the rate of primordial star formation 
prior to their assembly. In order for the first galaxies to remain metal-free, star forma- 
tion in minihaloes must be highly suppressed, most likely by H2-dissociating Lyman- 
Werner (LW) radiation. We show that the build-up of such a strong LW background 
is hindered by two effects. Firstly, the level of the LW background is self-regulated, 
being produced by the Population III (Pop III) star formation which it, in turn, sup- 
presses. Secondly, the high opacity to LW photons which is built up in the relic H II 
regions left by the first stars acts to diminish the global LW background. Accounting 
for a self-regulated LW background, we estimate a lower limit for the rate of Pop III 
star formation in minihaloes at z <; 15. Further, we simulate the formation of a 'first 
galaxy' with virial temperature Tvir ^ 10"* K and total mass ^ 10^ Mq at z ^ 10, 
and find that complete suppression of previous Pop III star formation is unlikely, with 
stars of ^ 100 Mq (Pop III.l) and ^ 10 Mg (Pop III.2) hkely forming. Finally, we 
discuss the implications of these results for the nature of the first galaxies, which may 
be observed by future missions such as the James Webb Space Telescope. 
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1 INTRODUCTION 

How did the first galaxies in the Universe form? In the hi- 
erarchical picture of structure formation, the first galaxies, 
with masses of ^ 10* M©, were built up from the mergers 
of smaller dark matter (DM) minihaloes, with virial tem- 
peratures ^ 10* K, in which the first Pop III star formation 
may have occurred (e.g. Bromm et al. 1999, 2002; Abel et 
al. 2002). The feedback effects from the first stars forming in 
minihaloes, likely having masses of the order of 100 M© (e.g. 
Tan & McKee 2004; Yoshida et al. 2006), may thus have es- 
tablished the properties of the gas from which the first galax- 
ies formed (see e.g. Ciardi & Ferrara 2005). During their 
brief lives of ~ 3 Myr, very massive Pop III stars emit co- 
pious UV radiation which can ionize the primordial gas and 
destroy H2 molecules (e.g. Schaerer 2002), and upon their 
deaths they may explode as powerful supernovae, expelling 
the first heavy elements into their surroundings (Mori et al. 
2002; Bromm et al. 2003; Kitayama et al. 2005; Greif et al. 
2007). 
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The nature of the first galaxies may largely be deter- 
mined by the metallicity of the gas from which they form, 
as the character of star formation is predicted to transition 
from a massive Pop III initial mass function (IMF) to a 
low-mass dominated IMF when the primordial gas has been 
enriched to a critical metallicity Zcrit (e.g. Bromm et al. 
2001a; Schneider et al. 2003; Santoro & ShuU 2006; Frebel 
et al. 2007; but see Jappsen et al. 2007). In turn, the metal- 
licity of the protogalactic gas is dependent on the preceding 
Pop III star formation that takes place in minihaloes. While 
the supernova explosions of Pop III stars formed during the 
assembly of the first galaxies could drive the metallicity of 
these systems to super-critical levels (see e.g. Greif et al. 
2007; Karlsson et al. 2008), Pop III star formation in mini- 
haloes could, in principle, be largely suppressed due to ra- 
diative feedback effects (e.g. Haiman et al. 2000; Mackey et 
al. 2003). 

Much recent work has been devoted to studying the 
effects of both ionizing and molecule-dissociating, Lyman- 
Werner (LW), radiation emitted from local star forming re- 
gions on Pop III star formation in minihaloes (Glover & 
Brand 2001; Shapiro et al. 2004; Alvarez et al. 2006; Abel 
et al. 2006; Susa & Umemura 2006; Ahn & Shapiro 2007; 
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Yoshida et al. 2007; Johnson et al. 2007; Whalen et al. 2007), 
with results generally suggesting that local intermittent LW 
feedback does little to delay star formation, while ionizing 
radiation only delays star formation in minihaloes with rel- 
atively low-density gas that may not efficiently form stars 
even in the absence of radiative effects (e.g. Ahn & Shapiro 
2007; Whalen et al. 2007) . Thus, it appears hkely that inter- 
mittent local radiative feedback effects may not be decisive 
in suppressing the formation of Pop III stars during the as- 
sembly of the first galaxies (see also Wise & Abel 2007a) . 

The properties of the first galaxies are, however, more 
likely related to the global LW background, as it has been 
shown that a strong, persistent LW background can sup- 
press star formation in the minihaloes which eventually are 
merged to form the first galaxies (e.g. Dekel & Rees 1987; 
Haiman et al. 1997; Bromm & Larson 2004). If star forma- 
tion in these small systems can be suppressed effectively by 
the LW background, the first protogalaxies are much more 
likely to be kept metal-free. However, if star formation is not 
easily suppressed in minihaloes, either because the LW back- 
ground is weak or if star formation can continue unimpeded 
despite strong LW feedback, then the first galaxies will likely 
form from already metal-enriched gas. Many studies have 
sought to estimate the global LW background during the 
epoch of the first stars (Haiman et al. 1997, 2000; Ricotti 
et al. 2002a,b; Glover & Brand 2003; Yoshida et al. 2003; 
Wise & Abel 2005; Greff & Bromm 2006). Additionally, the 
impact of a constant LW background on Pop HI star for- 
mation has been investigated, suggesting that, while global 
LW feedback does not completely suppress star formation 
in minihaloes, it can delay such star formation (Machacek 
et al. 2001; Mesinger et al. 2006; O'Shea & Norman 2008; 
Wise & Abel 2007b). 

In this paper, we investigate the impact of the build- 
up of a LW background on star formation in systems which 
later evolve into the first galaxies at redshift z J; 10. In Sec- 
tion 2, we develop analytical models which self-consistently 
couple the evolution of the LW background to that of the 
Pop HI star formation rate (SFR) aX z ^ 15. We describe 
our implementation of the LW background in cosmological 
simulations in Section 3. We report the outcome of simula- 
tions designed to test the hypothesis of a self-regulated LW 
background in Section 4, while in Section 5, we simulate the 
assembly of the first galaxy under the influence of such a 
self-regulated LW background. Finally, in Section 6 we dis- 
cuss the implications of our findings for the galaxy formation 
process, and we conclude in Section 7. 



2 THE LYMAN- WERNER BACKGROUND 

The LW background radiation field is produced by sources, 
either stars or miniquasars, at cosmological distances, be- 
ginning with the formation of the first stars at the end of 
the cosmic dark ages. This radiation field acts to delay star 
formation by dissociating H2 molecules, which are the pri- 
mary coolants allowing primordial gas to collapse and form 
stars in minihaloes. Thus, to fully address both the produc- 
tion of the LW background radiation and the Pop III star 
formation which is sensitively coupled to it would require 
a cosmological simulation resolving the collapse of gas into 
minihaloes over enormous volumes of the Universe. As car- 



rying out simulations over such large cosmological scales, 
while simultaneously resolving the collapse of gas in mini- 
haloes, is still prohibitively expensive, we revert to simple 
analytical estimates for the self-consistent build-up of the 
LW background which draw on the results of detailed nu- 
merical simulations of Pop HI star formation and radiative 
feedback. The simple self-consistent analytical model that 
we develop in this Section may serve as a first order ap- 
proximation to a more detailed self-consistent model for the 
build-up of the LW background and of the Pop III star for- 
mation rate in minihaloes at redshifts z J; 15. 

2.1 The Critical LW Flux 

We first estimate the LW fiux required to significantly delay 
the formation of Pop III stars in minihaloes. We note that 
for a typical case, following the extensive results of detailed 
simulations of Pop HI star formation (see e.g. Bromm et al. 
2002; Yoshida et al. 2003), the primordial gas first collapses 
adiabatically into a minihalo until it reaches a temperature 
T ^ 10'^ K and density n ^ \ cm^^. As the evolution of the 
gas is adiabatic up to this point, the cooling properties of 
the gas play a role in its evolution only at higher densities. 
In particular, for the gas to continue collapsing on the way 
to forming a Pop HI star, in accordance with the results 
of detailed simulations, a sufficient fraction of the coolant 
H2 is required (Tegmark et al. 1997). In the presence of 
a LW radiation field, however, H2 may be destroyed at a 
rate higher than it can be produced, effectively delaying the 
collapse of the gas. 

We may estimate the fiux of the radiation field required 
to achieve this delay by comparing the formation time of H2 
molecules, at T --^ 10^ K and n ~ 1 cm""^, to the dissociation 
time of H2 molecules as determined by the LW flux. For the 
dissociation time, we have (Abel et al. 1997) 

tdi.s ~ 3 X lO^Lw yr , (1) 

where Jlw is in units of 10^^^ erg s^^ cm^^ Hz~^ sr~^. 
Equating to the formation time found in our simulations. 
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8 X 10 yr, yields idis; 
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when the LW flux 



is JLw.crit ^ 0.04, which we deflne as the critical flux nec- 
essary to signiflcantly delay Pop III star formation in mini- 
haloes, roughly consistent with numerous simulation results 
(see Machacek et al. 2001; Yoshida et al. 2003; Mesinger et 
al. 2006; O'Shea & Norman 2008; Wise & Abel 2007b). 

2.2 The Maximum Background Flux 

As a first step in developing our estimate of a self-consistent 
Pop HI star formation rate and LW background flux, we cal- 
culate the maximum value for Jlw, as a function of redshift 
z, that we expect to be produced by Pop III stars form- 
ing in minihaloes at 2 ^ 15. We use the Sheth-Tormen for- 
malism (see Sheth et al. 2001), with the cosmological pa- 
rameters derived from the Wilkinson Microwave Anisotropy 
Probe ( WMAP) third year data (Spergel et al. 2007), to flnd 
the mass fraction collapsed in minihaloes at these high red- 
shifts. In comparison, the collapsed fraction in more massive 
haloes is very small, allowing us to neglect star formation 
that may be taking place in systems with virial tempera- 
tures Tvir ^ W K (see Section 2.3). Most studies of pri- 
mordial star formation in minihaloes find that single stars 
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Figure 1. The average star formation efficiency /* for Pop III 
star formation in minihaloes, assuming that a single 100 M© star 
forms in each collapsed minihalo with a virial temperature 2 X 10^ 
K < Tvir < 10* K. Here /* is defined as the fraction of collapsed 
baryons converted into Pop III stars. The average is carried out 
over the Sheth-Tormen mass function. 



with masses of the order of 100 M© form in these systems 
(Abel et al. 2002; Bromm & Loeb 2004; Yoshida et al. 2006; 
Gao et al. 2007). To estimate the maximum star formation 
rate we assume that each minihalo with a virial temperature 
Tvir ^ 2 X 10^ K (see e.g. Yoshida et al. 2003) hosts a 100 
Mq Pop III star which lives for t, — 3 Myr (e.g. Schaerer 
2002). The comoving star formation rate of Pop III stars in 
minihaloes, in the absence of the negative feedback due to 
LW radiation, can then be estimated to be 



SFRiii, 
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10- 
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16 
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where /, is the fraction of collapsed baryons converted into 
Pop III stars, averaged over the Sheth-Tormen mass func- 
tion, and \dF/dz\ is the change in the collapse fraction in 
minihaloes with redshift derived from the Sheth-Tormen for- 
malism. We calculate /» assuming that every minihalo with 
2 X 10^ K < Tvir < 10* K hosts a single Pop III star with a 
mass of 100 M0 , and the resulting evolution of /, with red- 
shift is shown in Fig. 1. Due to the increase of the average 
mass of such minihaloes with decreasing redshift, the frac- 
tion of collapsed baryons which goes into stars drops with 
redshift. For the following argument, it is convenient to com- 
pute a physical (proper) number density of Pop III stars in 
minihaloes as a function of redshift z, which we find to be 
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assuming a constant Pop III mass of 100 M0. 

To calculate the LW background generated in the lim- 
iting case that every minihalo hosts Pop III star formation, 
we next estimate the maximum physical distance traveled 
by a LW photon before being absorbed by a Lyman series 



transition in atomic hydrogen, which will be abundant at 
these epochs, when the Universe has not yet been substan- 
tially reionized (e.g. Haiman et al. 1997; Mackey et al. 2003). 
This distance corresponds to that of light traveling for a 
time needed to redshift the Ly/3 line of neutral hydrogen, at 
12.1 eV, to the low-energy end of the LW band at 11.2 eV. 
At this point essentially all LW photons, with energies be- 
tween 11.2 and 13.6 eV, will have been absorbed by neutral 
hydrogen. The distance to the LW horizon is thus 



10 



(^) 



Mpc 



(4) 



again given in physical units. We note that this distance is 
an upper limit to how far a LW photon can propagate; for a 
more detailed calculation of the frequency-dependent prop- 
agation of LW radiation see Haiman et al. (2000) . Using the 
number density of Pop III stars in equation (3) , we integrate 
the flux from all the stars within a sphere of radius rmax, ob- 
taining the following estimate for the maximum background 
LW flux: 
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( f* 
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(5) 



Following Johnson et al. (2007) , we here assume that the LW 
photons are emitted with a blackbody spectrum at 10^ K 
from a 100 Mq Pop III star, as modeled by Bromm et al. 
(2001b). In Fig. 2, we show the redshift evolution of JLW.max, 
along with the corresponding SFRni.max, assuming that /« 
evolves as shown in Fig. 1. 

We note that this idealized approach somewhat overes- 
timates the LW background from 100 M0 primordial stars, 
as many of the LW photons will be absorbed over shorter 
distances by the higher order Lyman series transitions corre- 
sponding to higher LW photon energies (e.g. Haiman et al. 
2000). Nevertheless, we expect that the above calculation 
gives a useful estimate for the maximum Jlw generated by 
Pop III star formation in minihaloes. 



2.3 Self-consistent Models of the LW Background 

In our calculation of the maximum LW background, we as- 
sumed that the Pop III SFR was independent of the level of 
the LW background generated. However, as shown in Fig. 2, 
the value of Jlw, max can exceed that necessary to signifi- 
cantly delay star formation (e.g. Machacek 2001; Mesinger 
2006; O'Shea & Norman 2008; Wise & Abel 2007b). There- 
fore, the model presented in Section 2.2 is inconsistent in 
that the SFR required to produce the maximum LW back- 
ground cannot be maintained under such a strong back- 
ground. Therefore, we now modify our model to more ac- 
curately couple the SFR to the LW background. 



2. 3. 1 Self -regulation 

The sources that produce the LW flux at redshifts z ^ 15 
may be stars or miniquasars, in principle. However, the ac- 
cretion onto massive black holes resulting from the collapse 
of Pop III stars is not efficient in general (Johnson & Bromm 
2007; Pelupessy et al. 2007), although it may be in relatively 
rare cases (see Li et al. 2007), and so the first miniquasars 
were likely not an important source of LW radiation in the 
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Figure 2. The Pop III star formation rates (top panel) and the 
corresponding LW baclcground fluxes {bottom panel) for the three 
cases discussed in Sections 2.2 and 2.3. The maximum possible 
LW background, JLW.maxi is generated for the case that every 
minihalo with a virial temperature T^ii > 2 X 10"^ K hosts a 
Pop III star, without the LW background in turn diminishing the 
SFR (dashed lines). The self-regulated model considers the cou- 
pling between the star formation rate, SFRni criti a^nd the LW 
background built-up by it, Jlw crit (solid lines). The minimum 
value for the LW background, Jlw shield: is produced for the case 
of a high opacity through the relic H II regions left by the flrst 
stars and a corresponding star formation rate SFRm shield: higher 
than SFRnijcrit and approaching SFRni,max (dotted lines). We 
expect that the true Pop III SFR in minihaloes has a value inter- 
mediate between SFRni,crit and SFRjjj shield- 



early Universe. We further assume that the dominant con- 
tribution to the LW background at z ^ 15 is due to Pop III 
star formation in minihaloes, and we here neglect atomic 
cooling haloes as additional sources. At these redshifts, ^ 
90 percent of the collapsed mass is in minihaloes which, if 
they host star formation, are expected to host single massive 
Pop III stars (Abel et al. 2002; Yoshida et al. 2006; Gao et al. 
2007). Therefore, if stars form in atomic-cooling haloes (e.g. 
Oh & Haiman 2002) with the same efficiency /« with which 
stars form in minihaloes, the vast majority of stars will form 
in minihaloes at z ^ 15, thus justifying our assumption. 

However, the star formation efficiency /, in atomic- 
cooling haloes is not well-constrained at z J; 15. Recently, 
Ricotti et al. (2008) have argued that /, can take a range 



of values, from ~ 10~ to ~ 10~ , roughly scaling as the 
square of the halo mass. Because most of the collapsed mass 
that goes into atomic-cooling haloes will reside in the least 
massive systems, we take this as evidence that the mass- 
weighted average of /« is likely near the lower end of this 
range, which is similar to minihalo efficiencies. In Fig. 3 we 
show the fraction of LW photons, /Lw.mini, emitted by Pop 
III stars formed in minihaloes with a star formation rate 
given by SFRm, crit, for various assumptions about the ef- 
ficiency of star formation and the IMF of stars formed in 
atomic-cooling haloes. While for minihalo star formation we 
assume that single 100 Mq stars form with an efficiency as 
given in Fig. 1, we consider values of/* = 10~^, 10~^, and 
10"'^ for atomic cooling haloes. Furthermore, we consider 
two IMFs for atomic-cooling haloes: a standard Salpeter 
IMF for the case of an already metal-enriched (Pop II) sys- 
tem, and a Pop III.2 IMF postulating that metal-free stars 
form with masses of the order of 10 Mq (Johnson & Bromm 
2006; see McKee & Tan 2007 for the terminology). The latter 
case is the predicted star formation mode in atomic cooling 
haloes, provided that they remain pristine (e.g. Johnson & 
Bromm 2006; Greif et al. 2008). For the Salpeter IMF, the 
number of LW photons emitted per baryon in stars is taken 
to be Tjuw = 4 X 10^, and for the Pop III.2 case, we assume 
7)LW = 2 X 10'' (see Greif & Bromm 2006). Because our 
model of self-regulated star formation assumes that mini- 
haloes produce the dominant portion of the LW background, 
it will only be valid for redshifts where /Lw.mini Jj 0.8. As- 
suming /, ^ 10^^ in atomic-cooling haloes, we thus expect 
that our self-regulated model of star formation in minihaloes 
is valid at redshifts 2 J; 15 for a Salpeter IMF, and at red- 
shifts 2; ^ 18 for a Pop III.2 IMF. As suggested by the results 
presented in Section 5, the IMF in such systems may transi- 
tion from a Pop III IMF to a Pop II IMF after the formation 
of only one or a few massive Pop III stars. 

With the assumption that Pop III stars in minihaloes 
contribute the dominant portion of the LW background, 
the production of the LW background becomes a purely 
self-regulated process: any increase in the star formation 
rate will translate into an increase in the LW background, 
which will in turn suppress the star formation rate, and vice 
versa. While a prohibitively expensive simulation would be 
required to self-consistently track the build-up of the LW 
background in this way (but see Yoshida et al. 2003), we 
suggest that a self-consistent flux would be similar to the 
value that we found for the critical LW flux, Jlw, crit. As ex- 
plained in Section 2.1, a LW flux much below this value will 
not efficiently delay the formation of molecules in minihaloes 
and so will only slightly suppress the star formation rate; in 
turn, a LW flux much above this value will delay the star 
formation rate substantially (e.g. O'Shea & Norman 2008; 
Wise & Abel 2007b). 

The star formation rate expected under the assump- 
tion that the LW background must not exceed Jlw, crit can 
be derived from the Sheth-Tormen formalism. We follow 
our treatment for the max;imum LW flux from Pop III star 
formation in minihaloes (Section 2.2), now requiring that 
Jlw, crit = 0.04 is not exceeded, as shown in Fig. 2. We then 
find the fraction /rog of minihaloes that must form stars in 
order to produce this self-regulated LW flux from the simple 
relation 
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Jrcg(Z) 



JlJW,cvit{z) 
JLW,max(2) 



(6) 



Assuming that Pop III stars form in minihaloes, if not pre- 
vented from doing so by negative feedback, as single, massive 
stars (e.g. O'Shea & Norman 2008; Wise & Abel 2007b), 
it is the number of minihaloes at a given redshift which 
determines the SFR. We thus infer that this critical rate, 
SFRiii_crit, is suppressed, at most, by a fraction /rcg com- 
pared to the maximum possible Pop III rate in minihaloes, 
SFRiii_niax. Because our value for JLw.max is an upper limit 
(see Section 2.2), the corresponding suppression factor, /rog, 
is a lower limit. 

The calculation presented here gives a simple estimate 
of the Pop III SFR at z J^ 15, accounting for the coupling of 
the SFR to the LW background generated by Pop III stars. 
In Section 4, we report on a simulation carried out to test 
the validity of this model as a first approximation to a more 
detailed self-consistent model of the Pop III star formation 
rate history. 



2.3.2 Shielding from Relic H II Regions 

While massive Pop III stars formed in minihaloes create 
large H II regions, with (proper) radii of a few kpc (Alvarez 
et al. 2006; Abel et al. 2006), these stars live for only Ss 3 
Myr, at which point the relic H II regions that are left be- 
hind begin to cool and recombine. The conditions in such 
relic H II regions are conducive to the rapid formation of H2 
and HD molecules (Nagakura & Omukai 2005; O'Shea et al. 
2005; Johnson & Bromm 2006; Yoshida et al. 2007). John- 
son et al. (2007) showed that the resulting opacity to LW 
photons, tlw, can be significant, exceeding unity across just 
a single relic H II region. Furthermore, these authors found 
that the high H2 fraction formed in these regions generally 
remains high, as molecules quickly reform in between the 
intermittent episodes of exposure to LW radiation produced 
by nearby, short-lived Pop III stars. 

To evaluate the effect of the opacity through relic H II 
regions on the build-up of the LW background, we now re- 
peat the calculation in Section 2.2, this time including a 
maximal optical depth to LW photons through relic H II 
regions, tlw, max. We employ an iterative procedure in order 
to determine self-consistently, at each redshift, the star for- 
mation rate, the optical depth to LW photons, and the LW 
background flux, all of which are dependent on each other. 
To begin, we initialize the calculation at z = 50, at which 
time we assume that tlw, max = 0, the star formation rate 
is equal to the maximal rate SFRm^max, and the LW back- 
ground flux is equal to the maximal value Jbw.max, and we 
integrate forward in time. 

To calculate the cosmological average tlw, max we adopt 
the results of Haiman et al. (2000), who find that at z = 15 
the average LW optical depth is tlw ~ 0.04 for a uniform H2 
fraction of /h2 = 2 x 10~®. Consistent with the results of 
Haiman et al., this cosmological average optical depth scales 
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where nH2 is the average number density of H2 in the inter- 
galactic medium (IGM). 

As Pop III star formation in minihaloes proceeds, the 



1.0 ™-^-...:^^T— 

0.8- 

0.6 

0.4- 

°-^ ~ Salpeter IMF 

0.0 ' I I I I I I I I I I 



1.0 







""""■■■-.. ~ '~ 


0.8 






1 0.6 
s 


- 




0.4 


- 




0.2 


- 


Pop I II. 2 IMF 


0.0 






4 





35 30 



25 



20 



15 



Figure 3. The fraction of the total LW flux contributed by 
Pop III stars formed in minihaloes /LW.mini in our self-regulated 
model (sec Section 2.3.1), for various choices of the star forma- 
tion efficiency /* in atomic-cooling haloes, which we assume are 
not subject to LW feedback. Atomic-cooling haloes in which Pop 
II stars form with a Salpeter initial mass function (IMF) (top 
panel) produce a smaller fraction of LW photons, for a given star 
formation efficiency, than do atomic-cooling haloes which host 
Pop III. 2 stars with masses of the order of 10 Mq (bottom panel). 
The solid lines correspond to a star formation efficiency of /* = 
IQ-"^ in atomic-cooling haloes, while the dashed and dotted lines 
correspond to /* = 10"^ and 0.1, respectively. Our self-regulated 
model assumes that Pop III stars formed in minihaloes emit the 
dominant portion of LW photons (/* > 0.8). For different choices 
of the IMF and the star formation efficiency, our model will be 
valid down to different redshifts. 



volume-filling fraction of relic H II regions increases with 
the number of Pop III stars which form and rapidly die. 
Johnson & Bromm (2007) found that the high residual free- 
electron fraction in relic H II regions can persist for of the 
order of 500 Myr, in turn driving the vigorous production of 
H2 molecules. These regions are thus expected to maintain a 
high H2 fraction for over a Hubble time, and their abundance 
is therefore assumed to be proportional to the cumulative 
star formation history and to only increase with time. We 
calculate the volume-filling fraction of relic H II regions, 
fun, by assuming that each Pop HI star which forms in a 
minihalo leaves behind a relic H II region. We find 
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where we have assumed that each Pop III star leaves behind 
a rehc H II region with a radius of ~ 4 kpc [(1 + z)/20\~^ in 
which the recombination time is longer than the lifetime of 
the star (see Abel et al. 2006; Alvarez et al. 2006; Johnson et 
al. 2007; Yoshida et al. 2007). The (physical) number density 
of relic H II regions, whii, taken to be proportional to the 
collapse fraction in minihaloes as explained above, is given 
by 



WhII ~ 10 /shield 



10-- 



F{z) 
6x 10- 



l + z 
16 



Mpc" 



(9) 



Here F{z) is the fraction of collapsed mass in minihaloes as 
a function of redshift, given by the Sheth-Tormen formal- 
ism, with the fiducial value of F{z) ~ 6 x lO-'' at 2 ~ 15, 
while F{z) ~ 10~^ at 2 ~ 20. We account for the suppression 
in the SFR due to a shielded LW background by including 
the factor /shield, defined as the ratio of the integrated SFR 
in the self-regulated model including an IGM optical depth, 
SFRiii^ahieid, to be estimated below, to the integrated max- 
imum SFR, SFRiii max, as described in Section 2.1: 



/shle 



/^°° SFRni,shieid(gOrf^' 
/f SFR„i,max(z')rf^' 



(10) 



We can then re- write the volume-filling factor of H II regions 
as 



/h 



' 0.1/shle 



/. 



5 X 10- 



F(z) 
6x 10- 



(11) 



For simplicity, we use a constant average value of /, = 5 x 
IQ-* in equation (11), consistent with the values shown in 
Fig. 1. 

To obtain the maximum average IGM optical depth to 
LW photons, we follow the results of Haiman et al. (2000) 
along with equation (7) , this time accounting for an elevated 
H2 abundance within a fraction /hii of the volume of the 
Universe. We thus find 



TLW.i 



= 4 



/hii 
0.5 



in, 
10- 



1 + z 
16 



(12) 



where we use /hj ~ 10^* as an upper limit for the average 
H2 fraction in relic H II regions and an average overdensity 
in these regions oi 5 = 10, as estimated from the simulation 
presented in Johnson et al. (2007). 

Next, we self-consistently calculate the minihalo Pop III 
SFR, again assuming that the LW background fiux must not 
exceed JLw.crit, but this time taking into account the optical 
depth TLw.max through relic H II regions. The fraction of 
minihaloes in which Pop III stars are able to form is then 
given by 



fr(z) 



JLW,crlt(z) 
J'LW,max(z) 



(13) 



analogous to the no-shielding case in equation (6). The intu- 
ition here is that to compensate for the loss due to the IGM 
opacity, a larger SFR is required to maintain a LW back- 



ground at the critical level, compared to SFRm,, 



how- 



ever, the SFR cannot exceed the maximum one. We thus 
find for the SFR of Pop III stars in minihaloes, including 
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Figure 4. The maximum IGM optical depth to LW photons, 
■TLW.max, due to the high fraction of H2 in the relic H II regions. 
This corresponds to SFRjji shield S'Ud a LW background flux of 
Jlw shield- Note that locally the optical depth can be greater than 
this average, exceeding unity through a single H II region, owing 
to self-shielding of H2 (see Johnson ct al. 2007). 



the effects of both a LW background and the IGM optical 
depth to LW photons 

SFRlll,shleld = /.(^)SFRlII.max . (14) 

Finally, the corresponding LW background flux is 



Jl 



= Jl 






,(15) 



as plotted in Fig. 2. The optical depth TLw.max produced in 
this model is shown in Fig. 4. We interpret Jlw, shield to be 
the minimal value for the LW background flux at z ^ 15. 
Correspondingly, we find that SFRm, shield is an upper limit 
to the Pop III SFR in minihaloes at these redshifts. 

In reality, the clustering of Pop III sources, coupled with 
the increasing SFR towards lower redshifts, will likely lead 
to the efficient photodissociation of H2 molecules in some 
relic H II regions even at z ^ 15, as well as to overlap 
between relic H II regions as /hii grows with time. Both 
effects will lower the opacity to LW photons through the 
IGM, and the global average LW background will eventu- 
ally exceed the level Jlw ~ 10""' needed for the dissociation 
of H2 molecules in these regions to outpace their formation 
(Johnson & Bromm 2007), although we note that we have 
not accounted for the local optical depth to LW photons 
through individual relic H II regions, which can exceed unity 
through a single relic H II region (Johnson et al. 2007). Fi- 
nally, we note that the emission of X-rays from miniquasars 
may, in principle, lead to a high opacity to LW photons, as it 
has been shown that the partial ionization from X-rays can 
enhance the generation of H2 in the IGM (e.g. Kuhlen & 
Madau 2005). However, as noted above, it is not clear that 
the accretion rates onto Pop III relic black holes are suffi- 
ciently high for this mechanism to dominate the generation 
of LW opacity (e.g. Pelupessy et al. 2007). 
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2.4 The Expected Pop III Star Formation Rate 

We expect the actual LW background at z ^ 15 to lie be- 
tween Jlw, shield and JLW.crit- The corresponding Pop III 
SFR in minihaloes which generates this LW background 
likely falls between SFRm^crit and SFRm, shield, since the 
shielding provided by relic H II regions can keep Jlw ^ 
JLW.crit- In turn, this will allow to keep SFRm ^ SFRm^crit, 
owing to the diminished strength of the negative LW feed- 
back. In the remainder of the present work we will take 
JLW.crit as a plausible maximum value for the LW back- 
ground flux, and we will use this value as an input in our sim- 
ulations of galaxy formation, thereby obtaining lower limits 
to the amount of star formation that may take place in the 
assembly of these systems at z ^ 15. We note that, although 
we follow the assembly of a galaxy to z ^ 12, the star for- 
mation that does take place in our simulation occurs at z ^ 
15. Thus, our simple assumption here that Jlw ~ JLW.crit 
even at z ;$ 15 does not impact the results of our simulations 
significantly. 



3 METHODOLOGY OF SIMULATIONS 

With the goal of estimating the degree to which Pop III 
star formation in minihaloes affects the properties of the gas 
from which the first galaxies form at z ^ 10, we have car- 
ried out two cosmological hydrodynamics simulations which 
track the detailed evolution of the primordial gas under the 
influence of a LW background radiation field. The first sim- 
ulation is designed to test the validity of our self-regulated 
model for the build-up of the LW background, and is de- 
scribed further in Section 4. The second sinmlation follows 
the assembly of a galaxy virializing at z ^ 10, assuming 
a critical LW background, JLw.crit. This simulation is de- 
scribed in Section 5. 

As with previous work, for our three-dimensional nu- 
merical simulations we employ the parallel version of GAD- 
GET (version 1), which includes a tree (hierarchical) gravity 
solver combined with the smoothed particle hydrodynamics 
(SPH) method for tracking the evolution of gas (Springel et 
al. 2001; Springel & Hernquist 2002). Along with H2, H+, H, 
II~, H"'", e~. He, He"*", and He+^, we have included the five 
deuterium species D, D"*", D~, HD and HD"*", using the same 
chemical network as in Johnson & Bromm (2006, 2007). The 
specific initial conditions for our simulations are described 
separately in Sections 4 and 5. 



3.1 Sink Particles 

The primary result obtained from the present simulations is 
the number of Pop III stars which are able to form in mini- 
haloes at z i^ 15 during the assembly of a galaxy virializing 
at z ^ 10. To track the number and location of the forma- 
tion sites for Pop III stars, we allow sink particles to form 
when the gas has collapsed to high densities. Specifically, 
sinks are created when the gas has collapsed to a threshold 
density rires, above which the mass resolution of our simu- 
lation is insufficient to reliably follow the continued collapse 
of the gas, as explained further in Johnson et al. (2007) . We 
note that the time at which a sink particle forms will pre- 
cede the time at which a star is able to form by roughly a 



free-fall time in gas with a density nres; while this can in- 
troduce an error of ^ 10 Myr in the precise timing of star 
formation, the conclusions that we derive from the results 
presented here are not significantly impacted. Once having 
formed, the sink particles act as tracers of star formation 
sites, each sink representing the formation site of a single 
Pop III star. We focus our attention on when and where 
these sinks emerge in order to estimate the number of stars 
formed, and to gauge the feedback effects from these stars. 

3.2 Implementation of Background Radiation 

Motivated by the considerations presented in Section 2 and 
further justified in Section 4, we have included in our sim- 
ulations the self-regulated model of the LW background ra- 
diation field, in the form of a global H2 dissociation rate. 
The LW background is taken to be a function of redshift, 
its evolution given by JLW,crit(z), as shown in Fig. 2. The 
dissociation rate is 



rCdiss (Zj 



10 JLW.crit 



(z)s- 



(16) 



In addition to the photodissociation of H2, the photode- 
struction of other species can, in principle, affect the cooling 
properties of the primordial gas. The main reaction produc- 
ing H2 being H~ -|- H ^ H2 -I- e~ (e.g. Peebles & Dicke 
1968; Lepp & ShuU 1984), the photodestruction of H~ is a 
process that might lead to suppression of the H2 fraction of 
the gas (Chuzhoy et al. 2007; see also Glover et al. 2006). 
However, we find the photodissociation of H^ to be unim- 
portant for the suppression of H2 in collapsing minihaloes. 
As discussed in Section 2.1, the gas in such minihaloes will 
adiabatically collapse to densities of ^ 1 cm~^ and temper- 
atures of ^ 10^ K, even if the fractions of H2 or H~ are low. 
Following this adiabatic collapse, the fate of the gas does 
depend sensitively on the abundances of these species, and 
it is only at this stage at which the photodissociation of H^ 
might begin to have an effect. We find that the chemical 
species at this stage, absent strong shocks or photoioniza- 
tion, are roughly in equilibrium when we include the effect 
of the LW background, and therefore suppression of the H2 
abundance can be directly tied to the abundance of H~ at 
this stage. At this temperature and density, we find that the 
H~ formation time is of the order of years. We find that 
the minimal photodissociation timescale of H~, produced 
for the case of the maximal LW background Jlw. max (see 
Section 2.2), is of the order of J; 10^ yr, roughly two orders 
of magnitude higher than the formation timescale. Thus, 
the H~ fraction will not be significantly affected by a back- 
ground radiation field, and in turn the abundance of H2 will 
not be significantly affected, allowing us to neglect the pho- 
todissociation of H~ as a factor in delaying star formation. 
Furthermore, we neglect the photodissociation of H^ due to 
radiation from local star formation for similar reasons, as 
explained in Johnson et al. (2007). 

We have also neglected Lya heating as a suppressant of 
minihalo star formation (see Chuzhoy & Shapiro 2007; Cia- 
rdi & Salvaterra 2007). For the Pop HI star formation rates, 
and corresponding production rates of UV photons, found in 
our simulations, Lya heating will not raise the temperature 
of the IGM appreciably, even for the maximal SFR consid- 
ered here, SFRm. max (see Section 2.2). Finally, we do not 
include a photoionizing background in our calculations. At 
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these high redshifts, the IGM is still predominantly neutral. 
Photoionization will thus generally occur only locally, when 
isolated Pop III stars form, an intermittent process which 
is not likely to dramatically affect the overall SFR (see e.g. 
Ahn & Shapiro 2007; Johnson et al. 2007; Whalen et al. 
2007; Wise & Abel 2007a; but see Ricotti et al. 2002a,b). 



4 TESTING THE SELF-REGULATED MODEL 

We have argued in Section 2.2 that the LW background and 
Pop III star formation rate at redshifts 2; ^ 15 are strongly 
coupled to one another, and that a self-consistent SFR is 
likely close to what is required to produce JLw.crit. To test 
this Ansatz, we have carried out a fiducial simulation with no 
LW feedback, and one with identical initial conditions, but 
with the critical LW background included. Comparing the 
resulting suppression of star formation in our simulations to 
the analytical ratio SFRiii.crit/SFRin^max, we can roughly 
gauge the validity of our simple self-regulated model. 



4.1 Initial Conditions 

We initialize these cosmological simulations according to the 
ACDM model at 2: = 100, adopting the cosmological pa- 
rameters rJm = 1 — f^A = 0.3, Q.B ~ 0.045, h = 0.7, and 
as = 0.9, close to the values measured by WMAP in its first 
year (Spergel et al. 2003). Here we use a periodic box with 
a comoving size L = 500 h^^ kpc, and particle numbers of 
A^'dm = A^SPH ~ 145"^ , where the SPH particle mass is mspH 
~ 710 M0. 

As we have already noted in Section 2, conducting a suf- 
ficiently large cosmological simulation to thoroughly follow 
the build-up of the LW background is prohibitively expen- 
sive. Our choice of initial conditions for the present simu- 
lation is thus purely a practical one, so that we may follow 
the evolution of a sufficient number of minihaloes in our 
cosmological box for our purpose, which is to compare the 
suppression of the rate of collapse of gas in minihaloes with 
the suppression of the star formation rate that is expected in 
our simple self-regulated model. Because this simulation is 
not the large-scale cosmological simulation that is necessary 
to fully test the self-regulated model, all we can derive from 
this simulation is whether or not the delay in the collapse of 
the few (~ 20) minihaloes that form in our small simulation 
box is roughly consistent with the global average suppres- 
sion of the star formation rate that we predict in the self- 
regulated model. Therefore, the results that we derive from 
this test hold no statistical significance which is dependent 
on the specific cosmological model that we use to initialize 
our simulation; we only seek to compare the isolated behav- 
ior of the minihaloes in our box with the general behavior 
we expect in the self-regulated model. A more detailed study 
of the self-consistent build-up of the LW background will be 
the subject of future work. 

The density above which resolution is lost in this sim- 
ulation is riros ~ 20 cm^*, as determined following Johnson 
et al. (2007). Therefore, sink particles are formed in mini- 
haloes which reach this density threshold, and we thereafter 
consider these minihaloes as hosting the formation of single 
Pop III stars. While we are not able to resolve the process 
of star formation itself, we believe this is a reliable indicator 
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Figure 5. The fraction of minihaloes which form stars when ex- 
posed to the critical LW background, JLW.criti normalized to the 
total number of stars formed absent a LW background, as a func- 
tion of redshift. Solid line: The fraction of Pop III stars expected 
to form in our analytical model for the critical star formation 
rate, SFRni,criti assuming the star formation efficiency /» shown 
in Fig. 1. Open triangles: The results from our simulation show 
the same general variation with z as the analytically-derived be- 
havior to within a factor of ^ 2, with the points from our simula- 
tion likely being overestimates due to our choice of WMAP first 
year initial conditions for our simulation. This agreement sup- 
ports the validity of our simple model for self-regulated Pop III 
star formation. 



of when star formation is possible, as gas collapsing in mini- 
haloes will generally only reach these densities when cooling 
is efficient and the collapse has ceased to be adiabatic (e.g. 
Omukai et al. 2005; Yoshida et al. 2006; Johnson et al. 2007). 



4.2 Results 

In Fig. 5, we compare our analytical prediction, described 
in Section 2.3.1, to our simulation results. The line in Fig. 5 
shows the analytically derived value for the fraction of stars, 
fcrit{z), which form by a redshift 2 in the self-regulated 
model out of the total number of minihaloes which could 
host star formation absent any negative LW feedback. We 
define this fraction as the ratio of the integrated SFR in the 
self-regulated model to the integrated SFR in the model for 
maximal star formation, described in Sections 2.3.1 and 2.2, 
respectively: 



/crit — 



/^°°SFRiii,crit(2')rf^' 

/^°° SFR„I,n,ax(2')d^' 



(17) 



In this analytical treatment, we have assumed the value 
shown in Fig. 1 for the fraction of collapsed baryons in stars 
/, , as a function of redshift. 

The quantity /crit is calculated as the number of sink 
particles that form by a redshift 2 in the simulation including 
the critical LW background, JLw.crit, divided by the number 
that form in the simulation without any LW background in- 
cluded, and is plotted in Fig. 5. As expected, likely owing 
to the relatively high bias in our simulations initialized with 
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WMAP first year parameters, the simulated values for /crit 
are higher than those from our analytical derivation using 
WMAP third year parameters. Nonetheless, the agreement 
between our simple analytical model and the simulations is 
broadly consistent, given the simplicity of our model. We 
thus conclude that, as a rough estimate, a self-consistent 
value for the LW background at redshifts z ^ 15 is likely 
close the critical value, JLw.crit- We emphasize, however, 
that the effect of the opacity to LW photons through relic 
H II regions will be to lower Jlw and, simultaneously, to al- 
low for a higher SFR than that in our simple self-regulated 
model shown in Fig. 2. 



5 GALAXY FORMATION IN THE EARLY 
UNIVERSE 

Assuming that the photodissociation of H2 molecules by a 
global LW background is the dominant factor in regulating 
Pop III star formation in minihaloes at redshifts z ^ 15, we 
now employ our simple model for a self-regulated LW back- 
ground to investigate the impact of Pop III star formation 
on the assembly of a galaxy at z ^ 10. We choose to define 
a galaxy as a star forming system hosted by a DM halo with 
a virial temperature of ^ 10* K, corresponding to a halo 
with total mass <; 10* M© at z <; 10 (e.g. Barkana & Loeb 
2001), as gas in such a system can cool via atomic hydrogen, 
can likely host continuous star formation, and can retain gas 
against being blown out by supernovae or by photoheating 
(e.g. Mori et al. 2002; Kitayama & Yoshida 2005; Read ct 
al. 2006). Here we report the results from our simulation of 
the formation of a galaxy which virializes at redshift z ~ 12, 
focusing on the modes and amount of star formation that 
may take place during its assembly. Additional aspects of 
the formation of such a galaxy are discussed in a companion 
paper (Greif et al. 2008). 
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Figure 6. The properties of the primordial gas just before the 
formation of the first sink particle at 2 = 16.3, as a function 
of distance from the sink: density (top-left), temperature (top- 
right), free electron fraction /^- (bottom-left), and H2 fraction 
/h2 (bottom-right). Note the enhanced free electron fraction in 
regions heated by the virialization of the host halo. Also, note that 
a second star later forms at 100 pc from the first in this halo, the 
second halo having collapsed to densities above n ~ 500 cm~'^ 
when the first star forms, densities which yield the evolution of 
the gas in the second halo relatively impervious to the radiative 
feedback from the first star (Whalen et al. 2007; Susa & Umemura 
2006; Ahn & Shapiro 2007). Finally, note that the H2 fraction is 
considerably lower than the values generally found for collapsing 
minihaloes in calculations not including a LW background. This 
leads to a higher temperature of the dense gas in the core of the 
minihalo compared to the canonical case absent a LW background 
(e.g. Bromm et al. 2002). This may lead to a larger mass for the 
Pop III star that forms from this gas. For reference, the resolution 
length in our simulation is ~ 3 pc (physical). 



5.1 Initial Conditions 

For our simulation of the assembly of a 'first galaxy', we have 
employed nmlti-grid initial conditions which offer higher res- 
olution in the region where the galaxy forms (e.g. Kawata & 
Gibson 2003; Li et al. 2007). We again initialize the simula- 
tion according to the ACDM model at z = 100, adopting the 
cosmological parameters Q.m = 1 — f^A = 0.3, Os = 0.045, 
h — 0.7, and erg = 0.9, close to the values measured by 
WMAP in its first year (Spergel et al. 2003). Here we use 
a periodic box with a comoving size L — 1 h~^ Mpc for 
the parent grid. This simulation uses a number of particles 
A'dm ~ NspYi ~ 1.05 X 10®, where the SPH particle mass is 
w-SPH ~ 120 M0 in the region with the highest resolution. 
For further details on the technique employed to generate 
our multi-grid initial conditions, see the companion paper 
by Greif et al. (2008). The density above which resolution is 
lost is riros ~ 10^ cm~^, as determined following Johnson et 
al. (2007). As described in Sections 3 and 4, sink particles 
are formed in minihaloes which reach this density threshold, 
and we thereafter consider these sinks to contain single Pop 
HI stars. 



5.2 Formation of the First Star 

The first star forms at redshift z ~ 16.3 inside the minihalo 
which is the most massive progenitor of the halo eventually 
hosting the galaxy. The mass of the DM halo hosting the 
first star is 7x 10® M©, while the virial temperature of the 
halo is 6 X 10^ K. This is comparable to the mass of the 
halo when star formation begins found in other studies for 
the same LW background and at similar redshift (O'Shea & 
Norman 2008). We note that this DM halo is approaching, 
but has not reached, virial equilibrium when the first star 
forms, having a gravitational potential energy of ~ 1.2 x 
10^^ erg and a kinetic energy of ~ 1.1 x IQ^^ erg. 

The properties of the primordial gas just before the star 
forms are different from the canonical case of Pop HI star 
formation in the absence of a LW background, as shown in 
Fig. 6. In the dense core of the halo, the H2 abundance is 
;$ 10^'', well below the canonical value of ;$ 10~^ at n ~ 
10"^ cm^"^, owing to the photodissociation of H2. In turn, 
the temperature of the gas in the core is ^ 500 K, markedly 
higher than in the case without photodissociation, for which 
temperatures are generally ~ 200 K. We note that the cen- 
tral temperature that we find in the protostellar core is sim- 
ilar to what is found in the simulations of O'Shea & Norman 
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(2008) for a comparable value of the LW background flux; 
however, these authors flnd a higher central H2 fraction of 
^ 10"*, likely because they follow the collapse of the gas to 
higher densities, n > 10* cm~"^, at which H2 formation may 
be faster. The higher temperatures of the collapsing gas will 
likely translate into higher accretion rates onto the proto- 
stellar core than in the case without photodissociation (e.g. 
Omukai & Palla 2001, 2003; Bromm & Loeb 2004; O'Shea & 
Norman 2008), and so possibly into a higher final mass for 
the star. If the mass of the star is pushed above ~ 260 Mq, 
or if it remains below ~ 140 M0 , it may collapse directly to 
form a black hole at the end of its life without yielding a su- 
pernova explosion, thereby keeping its surroundings largely 
metal-free (e.g. Fryer et al. 2001; Heger et al. 2003). 

The mass of the DM halo hosting this star is roughly an 
order of magnitude below that required for the virial tem- 
perature of the halo to exceed the atomic-cooling threshold 
of 10"^'^ K, and is likewise well below the mass that we ex- 
pect for the first galaxies. However, star formation is still 
able to take place at this relatively early time, despite the 
effect of the global LW background. We conclude that it is 
likely that at least a single Pop III star forms by z ~ 15 in 
a DM minihalo which is the most massive progenitor of an 
atomic-cooling halo hosting a galaxy at z ^ 10. The impact 
that such a star may have on the formation of the galaxy 
clearly depends sensitively on the final mass of this star. 

5.3 Local Stellar Feedback on Minihaloes 

We find that there is a second star-forming minihalo merely 
~ 100 pc from the site of the formation of the first star in 
our simulation, as can be seen in Fig. 6. Because of this prox- 
imity, we expect that any radiative, mechanical, or chemical 
feedback from the first star may dramatically impact the 
evolution of this neighbouring halo. Numerous recent stud- 
ies have discussed these feedback effects in detail (e.g. Ciardi 
& Ferrara 2005), and we here draw on them to understand 
how this stellar feedback may impact the formation of a 
galaxy in the early Universe. 

Regardless of whether the first star becomes a super- 
nova or collapses directly to form a black hole, the neigh- 
bouring halo will be subject to the ionizing and molecule- 
dissociating flux from this star (e.g. Shapiro et al. 2004; Susa 
& Umemura 2006; Ahn & Shapiro 2007; Johnson et al. 2007; 
Yoshida et al. 2007; Whalen et al. 2007). As the neighbour- 
ing halo will have already collapsed to densities ^ 500 cm""^ 
by the time the first star forms, the radiation from the first 
star is not likely to delay the collapse of the second halo 
significantly, as at these densities the core of the halo will 
not be photoevaporated and the molecules in the dense core 
may only be weakly dissociated (e.g. Susa 2007; Whalen 
et al. 2007). Therefore, it appears that a second star will 
form in the neighbouring minihalo within ^ 10 Myr after 
the formation of the first star (see also Wise & Abel 2007a). 
If the first star collapses directly to form a black hole, the 
formation of the second star will thus not be significantly 
influenced by the first star. 

The type of star that forms in the neighbouring halo 
is not entirely clear, however, if the first star becomes a 
supernova, owing to the mechanical and chemical feedback 
that this would entail. As shown in Greif et al. (2007), the 
velocity of the shock produced by a Pop III supernova at a 



distance of 100 pc from the progenitor is Vsh ^ 100 km s"^. A 
shock propagating at such velocities can strongly ionize the 
primordial gas (e.g. ShuU & McKee 1979), thereby allowing 
the formation of a high fraction of molecules and enhancing 
the cooling properties of the gas (e.g. Shapiro & Kang 1987; 
Mackey et al. 2003; Salvaterra et al. 2004; Machida et al. 
2005). In particular, high fractions of HD can be generated 
behind such strong shocks, allowing the post-shock gas to 
cool to the temperature of the cosmic microwave background 
(CMB) and perhaps to form metal-free stars with masses of 
the order of 10 M0, so-called Pop III. 2 stars. We note that 
for this type of Pop III. 2 star formation to take place the 
pre-shock gas must be at sufficiently high densities, unlike in 
the supernova simulation of Greif et al. (2007), in which this 
mode of star formation did not occur in the low-density IGM 
(see also Bromm et al. 2003). However, the neighbouring 
minihalo in the present case is likely to be sufficiently dense 
to allow for star formation to take place in the post-shock 
gas. We thus conclude that the halo neighbouring the first 
star may host a Pop HI. 2 star, in the case that the first star 
explodes as a supernova. 

The mode of star formation in the neighbouring mini- 
halo may change from Pop HI to Pop II, if metal-rich su- 
pernova ejecta from the first star can efficiently mix with 
the primordial gas in the dense core of the minihalo (e.g. 
Bromm et al. 2001; Santoro & ShuU 2006; Schneider et al. 
2006). We can estimate the efficiency of this mixing by ap- 
plying the criterion for the operation of Kelvin-Helmholtz 
(KH) instabilities given by Murray et al. (1993) (see also 
Gen & Riquelme 2008; Wyithe & Gen 2007); 

gDr 



27r«2 



^1 



(18) 



where g is the gravitational acceleration at the outer radius 
of the dense core of the neighbouring minihalo, and D the 
density ratio of gas in the halo compared to the dense shell. 
For the dense core of the neighbouring minihalo, we have a 
mass of the order of 10"" Mq and a radius of r ~ 30 pc. We 
find D ^ 500 and v^h Ss 100 km s"^, using values for the den- 
sity and velocity of the shocked gas from Greif et al. (2007) . 
The left-hand side of equation (18) is thus J; 10"^. Fol- 
lowing the results of Gen & Riquelme (2008), this suggests 
that the dense gas in the core of the neighbouring minihalo 
will remain largely pristine and stable for many dynami- 
cal times, while the outskirts of the halo may become more 
mixed. It therefore appears unlikely that the neighbouring 
halo will give rise to Pop II stars, but will instead proba- 
bly host Pop HI star formation. This will have implications 
for the nature of the galaxy into which these minihaloes are 
finally incorporated. 



5.4 Self-shielding and Molecule Formation 

While minihaloes are, in general, optically thin to LW pho- 
tons, the greater gas mass in a halo hosting the formation of 
a first galaxy can lead to a large column density of H2 and 
so to self-shielding of the molecules in the central regions 
of the halo (e.g. Draine & Bertoldi 1996). Taking this into 
account, we would like to estimate the impact of the LW 
background on the molecule fraction in the central regions 
of a first galaxy. Assuming that H2 formation takes place 
through the reaction sequence 
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Figure 7. The fraction of molecules in the primordial gas of 
the protogalaxy at 2 ~ 12.5, as functions of density: H2 fraction 
/h2 (left), and HD fraction /hd (right). The free electron frac- 
tion greatly exceeds the primordial value of f^- ~ 10"'* through 
coUisional ionization in virialization shocks, leading to enhanced 
fractions of both H2 and HD, despite the presence of the LW back- 
ground radiation field. Accounting for self-shielding, the molecule 
fraction remains elevated at high densities n l 10^ cm~''. In the 
core of the halo, the HD fraction is thus expected to reach values 
large enough to allow the formation of Pop in.2 stars (sec Grcif 
et al. 2008). 



e -I- H ^ H +hv 



(19) 
(20) 



we can estimate the formation rate of H2 as the gas collapses 
to high density. In our simulation, we find that the electron 
fraction drops due to recombinations approximately as 



/c 



10" 



nn 



102 cm- 



(21) 



where nu is the number density of hydrogen nuclei in the 
gas. Following Bromm et al. (2002), this results in a forma- 
tion rate of H2 at densities nn <; 10^ cm^"^ of 



V dt )i 



fcH-nn/o 



10" 



fcH-rijj 



(22) 



where fcn- 
tion (19). 

We can next estimate the photodissociation rate of H2, 
including the effect of self-shielding, following Bromm & 
Loeb (2003a). Accordingly, we estimate the column density 
of H2, using only local quantities, as Nb_^ ~ O.l/ns'^Hichar, 
where the local characteristic length is given by 



Lc 



/ 3XMb \ 
V A-KfrLYinii J 



1/3 



(23) 



Here Mb ~ 10^ Mq is the total baryonic mass in the halo, 
X = 0.76 is the mass fraction of hydrogen, and mn is the 
mass of the hydrogen nucleus. From our simulation of the 
formation of a first galaxy, which does not include the effects 
of self-shielding, we find a minimum average molecule frac- 
tion of /hj ~ 10~^ at densities nu ^ 10^ cm~^, as shown 
in Fig. 6, when Jlw = JLW.crit. The photodissociation rate 
will be reduced due to self-shielding by a fraction 



/sh 



iVi 



Ho 



10"cm- 



(24) 



when Nn^ > lO^'' cm^^ (Draine & Bertoldi 1996). We thus 
find a minimum photodissociation rate of 



4/H2 
dt 



3 X lO-^^nn'/Vu: 



1/4 



Jlw 



(25) 



Balancing the formation rate and dissociation rates for H2, 
given by equations (22) and (25), respectively, we can solve 
for the density at which the dissociation rate becomes higher 
than the formation rate, marking the point in the collapse 
of the gas where H2 dissociation begins to substantially sup- 
press the formation of H2. We find this density to be 



Wdi! 



10** cm 



-3 / Jhw \ 
VO.04/ 



/h2 
10-' 



-1/2 



(26) 



For the values of Jlw = 0.04 and /h2 ~ 10^^ in our simu- 
lation, we thus find that H2 photodissociation becomes im- 
portant only at very high densities. At densities below ndiss, 
the formation rate of H2 is higher than the photodissociation 
rate, owing to the higher abundance of free electrons, which 
catalyze the formation of H2 according to equations (19) 
and (20). As the density increases, the recombination of free 
electrons deprives the gas of this catalyst, and the forma- 
tion rate of H2 molecules thus drops, becoming comparable 
to the photodissociation rate when the gas has collapsed to 
a density ridiss- We note that photodissociation may not be- 
come important at all in the evolution of the gas, if ridias ^ 
10* cm-^, at which densities three-body reactions can effi- 
ciently produce H2 (see e.g. Glover 2005). Furthermore, the 
photodissociation of H2 is not likely to affect the evolution 
of the gas for cases in which the photodissociation timescale 
is longer than the free-fall timescale, as the gas may then 
cool and collapse to form a star before a significant portion 
of the H2 molecules are photodissociated. 

Along with the likely elevated fraction of H2, owing to 
ionization of the gas in the virialization of the galaxy. Fig. 
7 shows that a high fraction of HD is also generated. We 
expect that, due to the substantial self-shielding of the gas 
in the central regions of the assembling galaxy, the high HD 
fraction will persist in the gas as it collapses to high density, 
likely leading to the formation of Pop III. 2 stars (see Greif 
et al. 2008). 



6 IMPLICATIONS FOR THE ASSEMBLY OF 
THE FIRST GALAXIES 

In the preceding sections we have found that, owing in part 
to the self-regulation of the LW background, the Pop III star 
formation rate in minihaloes in the early Universe remains 
substantial at redshifts z ^ 15. In particular, we have found 
that the LW background is likely too weak to completely 
suppress star formation in the progenitor DM minihaloes 
that merge to from the first galaxies at 2: ^ 10. Here we 
discuss various implications of these results for the nature 
of the first galaxies. 



6.1 The Protogalactic Star Formation Rate 

While, in principle. Pop HI star formation can be suppressed 
in minihaloes by a global LW background, we have shown 
that with a simple model for the self-consistent build-up of 
the LW background Pop III star formation is likely to take 
place in the assembly of a galaxy at 2; J; 10 (see Fig. 8). 
Noting that the LW background that we used in our simu- 
lation did not take into account the effect of the opacity to 
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Figure 8. The hydrogen number density and temperature of the gas in the region of the forming galaxy at z ~ 12.5. The panels show the 
inner ~ 10.6 kpc (physical) of our cosmological box. The cluster of minihaloes harboring dense gas just left of the center in each panel is 
the site of the formation of the two Pop III stars which are able to form in our simulation including the effects of the self-regulated LW 
background. The remaining minihaloes are not able to form stars by this redshift, largely due to the photodissociation of the coolant H2. 
The main progenitor DM halo, which hosts the first star at z ^ 16, by z ^ 12.5 has accumulated a mass of 9 X 10^ M0 through mergers 
and accretion. Note that the gas in this halo has been heated to temperatures above 10^ K, leading to a high free electron fraction and 
high molecule fractions in the collapsing gas (see Fig. 7). The high HD fraction that is generated likely leads to the formation of Pop 
III. 2 stars in this system. 



LW photons through relic H II regions in the early Universe 
(see Section 2.3), the negative feedback from the global LW 
background may be even weaker than we found in our simu- 
lation. This, in turn, implies that the number of Pop III stars 
forming in protogalaxies may be higher than we find from 
our simulation. A total of two minihaloes collapsed to form 
Pop III stars in our simulation employing the self-regulated 
model for the build-up of the LW background, allowing us 
to conclude that of the order of a few Pop III stars are likely 
able to form during the assembly of a galaxy at high redshift. 
Additionally, we have found that conditions are prime 
for the formation of Pop III. 2 stars, owing to the high 
HD fraction that is generated in the virialization shocks in 
such systems. Because the Pop III. 2 star formation mode in 
atomic-cooling haloes has not been studied with simulations 
to the same extent as the 'classical', minihalo (Pop III.l) 
case, it is not yet possible to predict the resulting star for- 
mation rate and efficiency (but see Greif & Bromm 2006). 
However, we expect that, with 10^ M0 of gas collapsing into 
such a system. Pop HI. 2 stars may form in clusters, given 
that these stars are expected to have masses of the order of 
10 M0, compared to the 100 Mq predicted for Pop III.l stars 
formed in minihaloes. The stellar feedback is then expected 
to be less severe, possibly enabling cluster formation. 

6.2 Supernova Feedback and Metallicity 

Each of the two star-forming haloes in our simulation is lo- 
cated within ~ 100 pc of the center of the main progenitor of 



the protogalaxy, allowing each star to affect the inner region 
of the galaxy through either radiative output or by exploding 
as supernovae. We can estimate the degree of metal enrich- 
ment that may take place should each of these stars explode 
as pair-instability supernova (PISN), expected to be a com- 
mon fate for massive Pop III stars (e.g. Heger et al. 2003), 
as follows. For a 10* Mq DM halo, we expect ~ 10^ Mq of 
gas to collapse into the virialized system. With ~ 100 M0 
in metals being expelled into the medium from two PISNe, 
we thus expect to have a mean metallicity in this system of 
~ 10"^ Zq. 

Interestingly, this is above the value for the critical 
metallicity for Pop II star formation, Zcrit ~ 10"'^'^ Zq, 
as given by Bromm & Loeb (2003b) (see also Bromm et al. 
2001a; Frebel et al. 2007; Smith & Sigurdsson 2007). This 
metallicity estimate assumes very efficient mixing of the 
metals with the virialized primordial gas. In reality, there 
will likely be pockets of less enriched gas in the galaxy, 
suggesting that in the first galaxies there may be regions 
in which Pop II star formation is enabled through metal 
cooling, as well as other, more pristine regions in which a 
top-heavy IMF persists (see Jimenez & Haiman 2006). Con- 
sidering that the SFR may in fact be larger than what we 
found in our simulation, owing to weaker LW feedback, we 
expect that the mean metallicity in the first protogalax- 
ies may easily exceed the critical metallicity necessary for 
low-mass Pop II star formation. If the mixing of these met- 
als is efficient, the first galaxies may already host primar- 
ily Pop II star formation with a preponderance of low-mass 
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stars, provided that the mass of the stars formed in these 
systems is not dictated instead by other effects, e.g. the in- 
ability of the gas to cool below the temperature of the CMB, 
TcMB = 2.7 K (1 + z) (e.g. Larson 1998, 2005). 

However, the details of the transition from a top-heavy 
Pop III IMF to a low-mass Pop II IMF are not at present 
well understood, as it remains to understand the process 
by which the first metals, ejected into a hot expanding su- 
pernova remnant, are ultimately reincorporated into dense 
gas from which stars can form. The cooling time of the hot 
metal-rich ejecta from an energetic Pop III supernova is of 
the order of 10* yr (Greif et al. 2007). This is ample time 
for multiple Pop III stars to form from cold pockets of pri- 
mordial gas and explode as supernovae, allowing to enrich 
the gas with the yields from multiple supernovae before it 
re-collapses to form the first Pop II stars. Such a mixing 
of SN yields would dilute any pure PISN signature from a 
single massive Pop III star, and may explain the lack of 
observed metal-poor stars exhibiting a clear PISN signature 
(see Karlsson et al. 2008). An additional complication might 
arise from the presence of cosmic rays (CRs) produced by the 
first supernovae. Heating and ionization due to CRs could 
have an important efi'ect on the early star formation process 
(e.g. Jasche et al. 2007; Stacy & Bromm 2007), although 
these authors find that the main effect of cosmic rays is to 
enhance molecule formation in the central dense regions of 
minihaloes, which only strengthens our main conclusion that 
Pop HI star formation in minihaloes is not easily suppressed 
during the assembly of the first galaxies. 

6.3 Black Hole Formation 

A fraction of the first stars likely collapsed directly to form 
black holes, thereby locking up the majority of the met- 
als produced in their cores (see e.g. Madau & Rees 2001; 
Fryer et al. 2001; Heger et al. 2003). In a given protogalac- 
tic system, there is thus a probability that the primordial gas 
will not be enriched by supernovae even if stars are able to 
form within this system. Without knowing the IMF of Pop 
HI stars, however, we can only guess at what may be the 
fraction of Pop HI stars which collapse to form black holes 
without ejecting appreciable amounts of metals. Heger et al. 
(2003) give the mass range for Pop HI stars that collapse di- 
rectly to black holes as ~ 40 to ~ 140 M© and ^ 260 M©, 
while that over which a PISN takes place is given as ~ 140 
to ~ 260 Mq, for stellar models which neglect the effects of 
rotation. 

Given that Pop HI stars formed in minihaloes are char- 
acterized by a top-heavy IMF with typical masses of the 
order of 100 Mq, it seems reasonable to assume that the 
probability for a Pop HI star to collapse directly into a black 
hole, and thus locking up its metal-rich core, is of the order 
of 0.5. We note that this is consistent with the constraints on 
the fraction of PISNe from Pop HI stars reported by Karls- 
son et al. (2008). If we then take it that our simulation of the 
assembly of a ~ 10* M© galaxy is typical and use the fact 
that we find of the order of a few stars forming during this 
assembly process, we obtain a probability of the order of J^ 
0.1 that the first galaxy may be metal-free at the time of its 
formation. While this is only a rough estimate, the fact that 
we find little chance that star formation can be completely 
suppressed during the assembly of a galaxy at redshifts z ^ 



10 suggests that metal-free galaxies are probably a minority 
population of the first galaxies, with a majority of galaxies 
already hosting some amount of Pop II star formation. We 
emphasize, however, that Pop HI star formation may con- 
tinue to occur, even after the first supernovae, in regions of 
a galaxy which are not enriched with metals or in which the 
mixing of metal-rich supernova ejecta with primordial gas is 
inefficient. 

Finally, as discussed in Section 5.2, given that the LW 
background can suppress H2 cooling even in gas that does 
collapse to form Pop HI stars, the high accretion rates in the 
formation of these objects may allow for a large fraction of 
such Pop HI stars to have masses above ~ 260 Mq , allowing 
them to directly collapse into black holes and thus eject few 
metals. The LW background radiation may therefore act to 
keep galaxies metal-free even in cases where Pop HI star 
formation is not prevented. 



7 SUMMARY AND CONCLUSIONS 

We have investigated the impact of the H2-dissociating LW 
background radiation on the formation of Pop HI stars dur- 
ing the assembly of the first galaxies at z J; 12. To this end, 
we have constructed simple models of the build-up of this 
radiation field which are self-consistent, accounting for the 
coupling of the star formation rate to the generation of the 
LW background. We have motivated and, to the extent pos- 
sible, checked the self-consistency of a model in which the 
LW background does not exceed the critical value JLw.crit 
~ 0.04 at redshifts z ^ 15. Furthermore, we have elucidated 
the effect of a high opacity to LW photons through the first 
relic H II regions on the LW background, finding that the 
LW background flux may take values ^ JLW.crit at 15 ^ 2 
^ 18, if the cosmological average optical depth to LW pho- 
tons becomes sufficiently high. While our models are ideal- 
ized, they are informed by detailed cosmological simulations 
of Pop HI star formation and are consistent with previous 
results. We suggest that the actual LW background most 
likely assumes values intermediate between Jlw, shield and 
JLW.crit. However, detailed large scale simulations tracking 
star formation and accounting for the opacity of the IGM 
to LW photons will be required to better understand the 
true nature of the LW background at high redshifts. In ad- 
dition, observations of the cosmic near-infrared background 
could conceivably be used to place conservative upper lim- 
its on the level of the LW background at high redshifts, as 
LW photons would be reprocessed into Lya photons as they 
are redshifted with the expansion of the Universe, thus al- 
lowing them to contribute to the total Lya flux emitted at 
high z that may be observed today (e.g. Santos, Bromm & 
Kamionkowski 2002; Fernandez & Komatsu 2006). Taking 
into account the LW background, we find that the comov- 
ing rate of Pop HI star formation in minihaloes at 25 ^ z J; 
15 is likely <, 10"* M© yr"^ Mpc"^, with the upper limit 
corresponding to the case of a maximal opacity to LW pho- 
tons through relic H II regions at these redshifts. Overall, 
wo find that the global Pop HI SFR is likely decreased from 
the maximum SFR rate possible by a factor of, at most, ~ 3 
due to the LW background at 2 J; 15. 

We have simulated the assembly of a galaxy at z ^ 12, 
assuming the critical LW background in our self-regulated 
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model of global Pop III star formation. We find that of 
the order of a few Pop III stars are likely to form in the 
minihaloes that later merge to form the galaxy. Due to the 
chemical feedback that accompanies the explosions of Pop 
III stars as supernovae, possibly as PISNe, we conclude that 
most galaxies may already be enriched to a metallicity of 
the order of ^ 10^"' Zq when they are formed. If this is the 
case, then Pop II star formation may be widespread in these 
systems even at these early times. Alternatively, if the IMF 
of Pop III stars is such that a large fraction of them collapse 
directly to form black holes, thereby locking up the metals 
produced in their cores, then the interstellar medium of the 
first galaxies is more likely to remain metal-free. We find 
that a high fraction of the coolant HD is produced in the 
primordial gas during the virialization of the galaxy, despite 
the presence of the LW background, which likely leads to 
the formation of metal-free stars with masses of the order 
of 10 M© in these systems. If metal- free galaxies do exist, 
they may therefore be dominated by such a population of 
stars, perhaps the first stellar clusters (see also Clark et al. 
2008). More detailed simulations tracking the collapse and 
fragmentation of the gas, along with feedback from Pop III 
stars, will be necessary in order to test this scenario. 

The character of star formation clearly becomes more 
complicated in the first galaxies compared to the canonical 
picture of Pop III star formation in minihaloes, with the 
chemical feedback from the first stars, turbulence driven by 
rapid mergers and accretion of gas from the IGM (see Wise 
& Abel 2007c; Greif et al. 2008), HD cooling, CR heating 
and possibly magnetic fields all playing important roles. Fu- 
ture missions, such as the James Webb Space Telescope may 
thus be expected to observe some variety in the first galax- 
ies. The majority is likely already enriched with heavy ele- 
ments, hosting Pop II star formation, while a small fraction 
might still be largely metal-free, possibly hosting clusters 
of Pop III. 2 stars. The veil is about to be lifted, giving us 
access to the time when galaxy formation first began. 
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